Reduced density matrix of permutational invariant many-body systems 



o 
o 

(N 
> 

o 

OS 



43 ■ 

a,; 

^— > ■ 

c ■ 

=5 

cr 



> 



as 
o 



X 



Mario Salerno 1 and Vladislav Popkov 2 
1 Dipartimento di Fisica "E.R. Caianiello" and Istituto Nazionale di Fisica Nucleare (INFN), 
Gruppo Collegato di Salerno, Universitd di Salerno, 
via Ponte don Melillo, 1-84084 Fisciano (SA), Italy 
2 Dipartimento di Fisica "E.R. Caianiello" and Consorzio Nazionale per le Scienze Fisiche della Materia (CNISM), 
Universitd di Salerno, via Ponte don Melillo, 1-84084 Fisciano (SA), ttal^ 
(Dated: November 19, 2009) 

We consider density matrices which are sums of projectors on states spanning irreducible rep- 
resentations of the permutation group of L sites (eigenstates of permutational invariant quantum 
system with L sites) and construct reduced density matrix p n for blocks of size n < L by tracing 
out L — n sites, viewed as environment. Explicit analytic expressions of the elements of p n are 
given in the natural basis and the corresponding spectrum is derived. Results apply to all quantum 
many-body systems with permutational symmetry for which the mean field theory is exact. 



I. INTRODUCTION 

Reduced density matrices contain complete informa- 
tion about open quantum i.e. quantum systems in con- 
tact with an environment such as a thermal bath or a 
larger system of which the original system constitute a 
part (subsystem). In many cases one is interested in the 
spectrum of the reduced density matrix (RDM) because 
it reveals intrinsic properties (sometimes universal) of the 
subsystem. By definition, the spectrum of the RDM is 
real and nonnegative with all eigenvalues summing up to 
1. The relative importance of a state of the subsystem 
can then be measured by the weight the corresponding 
eigenvalue has in the RDM spectrum. Thus, for instance, 
the fact that the eigenvalues \ of the RDM for a one di- 
mensional quantum interacting subsystems decay expo- 
nentially with i implies that the properties of the subsys- 
tems are determined by only a few states. This property 
is crucial for the success of the density-matrix renormal- 
ization group (DMRG) method [l[ in one dimension. In 
two dimensions this property is lost 0] and the DMRG 
method fails. 

For a subsystem consisting of n sites (or n q-bits) the 
RDM is of rank 2" so that for large n the calculation of 
the spectrum becomes a problem of exponential difficulty. 
While the spectrum of the full RDM for subsystems with 
a small number of sites (e.g. n < 6) has been calculated 
the full RDM for arbitrary n, to our knowledge, is 
exactly known only for the very special case of systems 
of free fermions (see e.g. [j|). 

The aim of the present paper is to analytically calcu- 
late the elements of the RDM of permutationally invari- 
ant quantum systems of arbitrary size L, for arbitrary 
permutational symmetry of the state of the system (la- 
beled by an integer number < r < L/2) and for arbi- 
trary sizes n (number of q-bits) of the subsystem. We re- 
mark that the invariance under the permutational group 
physically implies that the interactions among sites have 



infinite range. From this point of view our results may 
apply to all quantum systems for which mean field the- 
ory becomes exact. As an example we consider a system 
of Heisenberg spins 1/2 on a full graph consisting of L 
sites, with fixed value of magnetization S z = L/2 — N. 
For this system we calculate the RDM for a subsystem 
of arbitrary n > 1 sites for arbitrary L, N. 

The plan of the paper is the following. The formula- 
tion of the main problem and the basic definitions are 
given Section 2. In Section 3 we discuss general proper- 
ties of the elements of the RDM while in Section 4 we 
use symmetry properties of the system to decompose the 
RDM into block diagonal form. The main results of the 
paper are presented in Sect. 5 in the form of a theorem 
giving the analytical expressions of the RDM elements 
for arbitrary L,n,r. For simplicity, we provide a proof 
of this theorem only in the thermodynamic limit. In Sec. 
6 we characterize the spectrum of the RDM and discuss 
some of its main properties. Finally, in the last section 
we briefly summarize the main results of the paper. 



II. MODEL EQUATION AND MAIN 
DEFINITIONS 



Consider a permutationally invariant system of L spins 
1/2 on a complete graph with fixed total magnetization 
S z = L/2 — N and described by the Hamiltonian 
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Here S = (S x ,S v ,S t ), S a = lE^of, with of Pauli 

L 

matrices acting on the factorized J\ space. This 

1 

Hamiltonian is invariant under the action of the symmet- 
ric group S l and conserves the total spin polarization S z , 
[H, S z ] = 0. A complete set of eigenstates of H are states 
\^L,N,r) associated to filled Young Tableau (YT) of type 



2 



{L~r,r} {N) (see 0), 

H\V L ,N,r) = E L ,N,r\^L,N,r), 



(2) 

E l ,n,t = \ (L - r + 1) + h{L - 2N) Y (3) 



.S':hI'L..N.,-) = ( ~-N j 



Here N = 0, 1, L determines possible values of the spin 
polarization and r takes values r = 0, 1, max(iV, L—N) 
(the explicit form of the state \^L,N,r) is given below in 
Eq ([25|) ). The degeneracies of the eigenvalues El,n,t are 
given by the dimension of the corresponding YTs, 



deg L , 



L 



L 

r - 1 



(5) 



Definition 1 Consider a set of vectors \*f? u ), u — 
l,...deg ir , forming an orthonormal basis in the 
eigenspace of H with eigenvalue pL,JV,r- We define the 
density matrix of the whole system as 



1 



c Sl, 



CTL,N,r 



deg z 



(6) 



&L,N,r possess the following properties: 
i) The matrix OL.jv.r has eigenvalues Ai = A 2 = ... = 
^de S L,r = (degL.r)" 1 - witn remaining 2 L - deg L r eigen- 
values all equal to zero. This follows from the fact that 
each vector \^ u ) is an eigenvector of ul.n.t with eigen- 



value 



Since the spectrum of ctl N r is real and 



nonnegative with all eigenvalues summing up to 1, the 
remaining 2 L — deg L r eigenvalues must vanish. 



ii) Matrix gl,n,t satisfies: ((JL,N,rY 



dcg £ 



This 



follows from the definition ([6]) and the orthonormality 
condition (* w |* u ) = S uw . 

iii) Introduce the operator Py, permuting subspaces i 

L 

and j of the Hilbert space Y[ §§0% on which the matrix 

l 

&L.N.r acts. Then [a, Pij] = for any i,j. 
Proof. 

Let us consider 



dG Sl,,r 



= — - — V I*' )(*' I 



(7) 



The vectors l^' 
Indeed, I*') 



because P„^ = Py, and (P;^ 2 



^d0gr 



y 



form an orthonormal basis. 

I. Now, the sum 



J2u=i ' r I I = Idc gL ,. is a unity operator in a sub- 

space of dimension deg L r , and therefore it does not de- 
pend on the choice of the basis. Note that vector \^' u ) be- 
longs to the same subspace as \^ u ), because permutation 



Pii only results in different enumeration. Consequently, 



Pij pPij 



or [p, P^} = 0. 



(8) 



The latter property implies that in Eq. ([6]) the sum over 
the orthogonalized set of basis vector in © can be re- 
placed by the symmetrization of the density matrix di- 
( 4 ) rectly, namely o L , N , r = 4 Y, P |*12...£>(*12...z|, where 



the sum is over all L\ permutations of indexes 1,2, ...L, 
and \^i2...l) is some unit eigenvector of H with eigen- 
value El,n,t- In particular, it is convenient to choose 

1*12.. .h) = Y^L.N.r), 



&L,N,r 



-y 

p 



I* 



L,N 



r)(* 



L.N.', 



(9) 



It is evident that such a sum is invariant with respect 
to permutations and that <Jl,n,t is properly normalized: 
Tra L ,N,r = 1. 

Definition 2 The Reduced Density Matrix (RDM) of 
a subsystem of n sites is defined by tracing out L — n 
degrees of freedom from the density matrix of the whole 
system: 



P(n) — Tr L ^ n UL. 



N,: 



(10) 



Due to the properties JHJ) and (fT0|) . p^ does not depend 
on the particular choice of the n sites, and satisfies the 
property {§} in its subspace (we omit the explicit depen- 
dence of pi n ) on L, N, r for brevity of notations). 

III. PROPERTIES OF THE RDM 

The RDM can be calculated in the natural basis by us- 
ing its definition in terms of observables: (/) = Tr(pr n \f) 
where / is a physical operator acting on the Hilbert space 
of the 2™ x 2™ subsystem. The knowledge of the full set 
of observables determines the RDM uniquely. Indeed, if 
we introduce the natural basis in the Hilbert space of the 

n 

subsystem, Yl <H)C2) the elements of the RDM in this 

fe=i 

basis are 

Jii ? -in _ /giii ? .~i„ \ = j, I gnis-tn \ (11) 



n« 

fe=i 



De*-* and a 2 x 2 matrix with 



elements (ej) = SikSji. The matrix eV*?"'*^ has only 
one nonzero element, equal to 1, at the crossing of the 
row 2 n - 1 i 1 + 2 n ~ 2 i 2 + ■■■ + i n + l and the column 2"~ 1 ji + 
2™~ 2 j2 + ••• + in + 1 (all indices i,j take binary values 
i k = 0,1 and j k = 0,1). To determine all the RDM 
elements one must find a complete set of observables and 
compute the averages (e£*? '"j™ )■ Note that a generic 
property of the RDM elements, which follows directly 
from ([8]) , is that any permutation between pairs of indices 
(im,jm) and (ik,jk) does not change its value, e.g. 



■3n 



= P 1211 - 



323l---3n 



P 



3-n.3\--3l 



= p 



3n3n-X—3\ ' 



(12) 
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Another property of the RDM follows from the S z invari- 
ance 



lll2...l n 



0, if ii + i 2 + ••• + i n 7^ h + h + ■•• + in- 

(13) 

Thus, for instance, the RDM for n = 2 has only three 
different nonzero elements, PooiPoi = Pio an( ^ Pu- It i s 
convenient to introduce the operators 







e =\ 1Q =* , e i = I o Of =a ' 



1 



1 




P , e x 




1 



= ft . 



If we represent a site spin up with the vector ( Q ) and a site 

spin down with the vector (°) then pk and ftfc are spin up 
and spin down number operators on site fc, while o~ , <r + 
represent spin lowering and rising operators, respectively. 
Thus, for instance, the observable {p\P2hj l hn...h n ) — 
PooiL'"! gi ves the probability to find spins down at sites 
3,4, ...n, and spins up at sites 1,2, while the observable 
( a t a 2 a 3 a 4 h5---h n ) = Piiooi"'! gi yes the probability to 
find spins down at sites 5,6, ...n, spin lowering at sites 
3,4 and spin rising at sites 1,2. Note that the latter 
operator conserves the total spin polarization since the 
number of lowering and rising operators is the same. Also 
note that the correlation functions with a non conserved 
polarization vanish, e.g. 



in accordance with (TT 



1)0111.. .1 
F11011...1 



0. 



IV. BLOCK DIAGONAL FORM OF THE RDM 

One can take advantage of the S z invariance (e.g. 
Eq. [T3")) to block diagonalize the RDM into independent 
blocks Bk of fixed polarization k = i\ + i% + ... + i n — 
ii+ia+-"+in (here k gives the number of spin up present 
in the subsystem). Note (see FigLU that the block di- 
agonal form in the natural basis becomes evident after a 
number of permutations of rows and columns of the RDM 
have been performed. Also notice that the n + 1 blocks 
correspond to the values s z = (n — 2k)/2, k = 0, 1, ...,n 
the subsystem polarization can assume. The dimension 
of a block Bk coincides with the number of possible con- 
figurations that k spin up can assume on n sites, e.g. 
dim Bk = (?) . One can easily check that the sum of the 
dimensions of all blocks gives the full RDM dimension, 
i.e. E fc dimS fc = 2". 



Blocks Bi. consist of elements e!- 



of 



the original matrix, with — ^2 j P — k and i p = 0, 1, 

l l 

j„ — 0,1. In its turn, all elements el 1 ®el 2 ®e\ 3 ®...®e ln 

*^ Jl .72 J3 J ti 

of the block Bk can be further block diagonalized accord- 
ing to the irreducible representations of the permutation 




FIG. 1: Left panel. Reduced Density Matrix pr 6 ) in the 
natural basis. Parameters are L = 12, N = 6, r = 3. 
Black boxes denote non zero elements. Right panel. The 
same matrix of the left panel after the following chain of 
permutations of rows and columns have been applied to it: 

f 22, 25 f 21, 34^20, 35 fl6, 37^15,41 P14, 33 fl2,21 ^8,49 Pr, 33 Pe, 17 Pi, 9- 

(Pi,j denote the operator which exchange first columns i 
and j and then rows i and j). Black boxes denote non zero 
elements. Block diagonal structure associated with values of 
k = 0, 1, 6 is evident. The single element present in blocks 
k = 0,6, is 1/924. Elements values inside other k-blocks are 
given in Fig. 



group S n of the subsystem (see also section 6 below) . In 
the natural basis, this block diagonalization is done ac- 
cording to the number of pairs of type (e^OeJ) present in 
the elements. In the following we denote by Gz the part 
of the block associated to elements with Z pairs (e\ £S> ej) 
in it. The sub-block Go of the block Bk is formed by 
the elements containing k terms e\ and (n — k) terms 
e§ in the product, i.e. e\ ® ...e\ ®e°® ... <g> e§ and all 
permutations. All such elements lie on the diagonal, and 
vice versa, each diagonal element of Bk belongs to Go- 
Consequently, the sub-block Go consists of (?) elements. 

The number of elements, degGi(fc), in the sub-block 
Gi is equal to the number of elements of the type e\®e\® 
e-J®e^(g)...<g)e^, such that l+0+h+i 2 + ■■■+i n -2 = k. 
Using elementary combinatorics we obtain: 



n-2 
k-1 



(14) 



Analogous calculations for arbitrary sub-block Gz yields 
r 2Z\ ( n\ fn-2Z^ 



deg G z (k) 



Z I \ 2Z \k-Z 



(15) 



From the restriction i p — j p — k we de- 
duce that the block Bk contains non-empty parts 
Go,Gi, G m i„(fc -n _fc), leading to the following decom- 
position: 



(16) 



min(A;,n — k) 

B k = |J G z . 

z=o 



Indeed, the normalization condition following from (|16p . 
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FIG. 2: Blocks fe = 1, 5 (left), k = 2,4 (center) and k = 3 
(right) appearing in the block diagonal form of p( 6 j given in 
the right panel of Fig[T] White boxes in the left panel de- 
note the element 1/154 while black boxes denote the element 
1/924. In the center panel white boxes denote the element 
5/1848 while black boxes, except for the diagonal elements 
which are equal to 5/308, denote the element —1/924. In the 
right panel white boxes denote the element 5/1386 except for 
the ones along the anti-diagonal which are all zero, while the 
black boxes denote the element —1/693, except for the ones 
along the diagonal which are all 5/231. 



gives 



min(/c,n— k) 

E 

z=o 



degG z (fc) 



(17) 



It is important to note that all elements of Gz are 
equal. This is a direct consequence of the property in Eq. 
([T!?)) . A graphical representation of the Bk block for a 
particular choice of L,N, r, is given in Fig|2] 



V. ANALYTICAL EXPRESSION FOR RDM 
ELEMENTS 

The main result of the paper is provided by the follow- 
ing theorem. 

Theorem Elements of RDM gz, belonging to a sub- 
block Gz of a block Bk of the reduced density matrix 
(1 1 0[) are given, for arbitrary L,N,r, n, by 



gz 



(L-n s 
UV-fc. 



Era=0 ( 



1 



-\ iL — N—r\ I r \ 
%) \ Z-m ) \m) 



(n) 



L-N\ 
Z ) 



(18) 



This expression defines completely all elements of the 
RDM in the natural basis. In practice, to find the el- 
ement P,Q of the RDM, {p( n ))pQ> we have to take the 
binary representation of numbers P — 1 and Q — l (which 
gives us the sets {i p } and {j p } respectively), find the 
corresponding number Z, and use (fT5|) . 

The proof of this theorem will be given in the ther- 
modynamic limit, L — > oo, ^ = p, £ = p,, in which the 



expression (jl8f simplifies to 

gz=P n - k (l- P ) k r ] Z , 



with 



(p-^)(i -p-p) 

p(l-p) 



(19) 



(20) 



We shall prove Eq. (pT9|) first for the case r = cor- 
responding to the symmetric ground states in (|21|) and 
then for the case r > corresponding to excited states. 



A. Case r = 0. 

If r = 0, the Young Tableau Yi,,iv> is nondegener- 
ate and the state of the initial system is pure: p = 
I* l,n){^ l,n\ with \^l,n) the symmetric state 



I* 



L,N) — 



-1/2 



(21) 



where the sum is over all possible permutations, and the 
prefactor takes care of normalization. Since all L sites are 
equivalent due to permutational invariance, any choice of 
n sites ii,l2, ■■■,i n within L sites gives the same RDM, 

It has been shown 



which we denote by p^^ : ° = Tr^- 



[Tfjj that p^'° takes form 



n f L — n\ 
L,N,0 ST^ \N-k) 



P(n) 



k=o (jv) 



|*n,fc)(*n,fe| 



In the natural basis the matrix elements of RDM are 
given by the above discussed values of observables. Using 
(|2"Tj) . one explicitly computes all RDM elements as 



( l,n,o\ ( l,n,o\ 

(/V) J =^(n) J 



_( L,N,0\ %J "' m _ x i+j + .. + .m \N-w) 



V+j'...+m' I L\ 

\n) 



(22) 

with w — i + j + ... + m (the sets ij...m and i'j'...m' are 
binary reprentation of numbers P — 1, Q — 1). So, all the 
elements of a given block Bk are equal. This property 
does not hold in the general case r > 0. Note that 
in the thermodynamic limit r\ = 1 , and / ( m) ~ * 

p n ~ k {l-p) k . 



B. General case r > 0. 

To calculate the RDM, we shall use the representation 
([9]) for the density matrix of the whole system a, rewrit- 
ing it in the following form 



P(n) 



Tr L - 



Tr L - 




n! (L-n)\(t) 



Pr„i Pi, 



\^L,N,r)(^L,N,r\ 



(23) 
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FIG. 3: Graphical representation of a filled Young Tableau 
{L — r, r}jv with L = 11, r = A,N = 6. The 1 inside a box 
denote a spin up (particle) while the denote a spin down 
(hole). 



Note that the L\ permutations can be done in three 
steps: first, choose at random n sites i\ ^ i% =/= ■■■ =/= i n 
among the L sites. There are (^) such choices. Then, 
permute the chosen n sites, the total number of such 
permutations being n\. Finally, permute the remaining 
L — n sites, the total number of such permutations being 
(L — n) ! . The latter step (c) under the trace operation is 
irrelevant because these degrees of freedom will be traced 
out. The operation permuting n sites commutes with the 
trace operation since Tr^_„ does not touch the respective 
subset of n sites. Consequently, (|23p can be rewritten as 



P{n) 



iH Tr ^«7lA \^L,N,r){^L,N,r\ 



(24) 

The filled Young tableau {L — r,r}jy consists of the 
antisymmetric part containing 2r sites and r spin up, 
and the symmetric part containing the remaining L — 2r 
sites and N — r spin up, see Fig. [3] (in the following we 
adopt an equivalent terminology appropriate for non spin 
systems by identifying spins up with particles and spins 
down with holes). The wave function \^l,n,t}] corre- 
sponding to this Young tableau factorizes into symmetric 
and antisymmetric parts 

* = \(t>12)®\foi)®-~®\4>2r-l,2r)®\^L-2r.N-r)2r+1.2r+2...L 

(25) 

with the antisymmetric part consisting of the first r fac- 
tors of the type 



|^ 12 ) = ^=(|10)l2-|01)l 2 ) 



(26) 



and with the symmetric part, \^ L-2r.N-r) , given by (|21[) . 
A general property of factorized states implies that if the 
global wave function is factorized, |<f>) = \ip)i\<f))ii and 
out of n sites of the subsystem, n\ sites belong to subset 
J, and the remaining n 2 — n — n\ sites belong to the 
subset II, then the reduced density matrix factorizes as 
well: 

P(n) =P(n 1 ) ( 27 ) 

To do the averaging, we note that among total number 
of choices ( L ) there are (a) ( L ~ 2r ) possibilities to choose 
n sites inside the symmetric part of the tableau, contain- 
ing N — r particles, (b) 2r(^~ 2 1 r ) possibilities to choose 
n — 1 sites inside the symmetric part of the tableau and 



one site in the antisymmetric part (c) ( 2 2 r ) {^~_ 2r ) possi- 
bilities to choose n — 2 sites inside the symmetric part 
of the tableau and two sites in the antisymmetric part 
and so on. The contributions given by (a) and (b) to the 
right hand side of for ' r are given, according to 
USB, by 



F \ F,M,0 
P(n) 



2r 



F 



1 



F,M,0 
P(n-l) 



(28) 



with F = L - 2r, M = N -r and with pi = 1/2 the 
density matrix corresponding to a single site in the an- 
tisymmetric part of the tableau. Brackets (.) — J^p 
denote the average with respect to permutations of n ele- 
ments. The contribution due to (c) to (|24|) splits into two 
parts since the ( 2 2 r ) possibilities to choose two sites in the 
antisymmetric part of the tableau consist of 4Q choices 
with two sites into different columns and the remaining 
r choices with both sites belonging to a same column. 
For the former choice, the corresponding density matrix 
pi , while for the latter case is given by 

, with 



f.m,o 
is P ( „_ 2) 65 pi 

F,M,0 „ 
P(n—2) ^ Pdsymnt 



Pa 



SJ/mm = |-L(10-01))(-^(10-01)| 



/0 



0\ 
1-10 
0-110 

Vo o oo/ 

(29) 

Proceeding in the same manner for arbitrary partitions 
of Z sites in the antisymmetric part of the tableau and 

n 



Z sites in the symmetric part, we get 

F 



L,N,r I L 
P(n) 



II ® Pas 



I min(2r,n) 



,Z-2i 



[Z/2] 
F,M,0 \ " 
P(n-Z) 2^ 




From this the general scheme for the decomposition of 
the general RDM becomes evident. In the above formula, 
Q 

the products Y[ with Q < i are discarded. The matrix 

i 

elements p^ 1 ' are given by (|22|) . 

For simplicity of presentation, we prove Eq. (|19p for 
the case Z = k and then outline the proof for arbitrary 
Z. 

In the thermodynamic limit one can neglect the differ- 
ence between factors like 4(!J) and ( 2 2 r ) in Eq. ([50)1 . The 
latter then can be then rewritten in a simpler form as 



L\ L.N,; 
P(n) 



F\ F,M.O 
n )P(n) 



F 
n — 1 



F.M.O 
P(n-1) 



3 PL 

(31) 



F 



F,M,0 
P(n-2) ' 



' Pi ® P I 



Note that one can omit all terms in (|50"|) containing 
pasymm since the respective coefficients correspond to 



6 



probabilities of finding two adjacent sites in the asym- 
metric part of the YT (proportional to r), which vanish 
in the thermodynamic limit, respect to the total num- 
ber of choices which is of order of r 2 . A sub-block Gz 
of a block k consists of all elements of the matrix 
having Z pairs of e\,e\ in its tensor representation, like 
e.g. (e$ <g> ej)® z ® e\ <g> eg <g> ... ® e]"~^, such that 
Z + ii + i2 + ■■■ + i n -2Z = k. The total number of elements 
gz C Gz in Pr^^ 1S equal to the number of distributions 
of Z objects e\, Z objects ej, and (k — Z) objects e\ on 
n places, given by 



deg G z = 



Z\Z\{k - Z)l(n-k - Z)\ 



(32) 



(this is another way of writing (|15|) V Each term W in 
the sum (|31[) after averaging will acquire the factor 



T(W) = 



deg G Z {W) 
degGz 



(33) 



where deg Gz(W) is a total number of gz elements in 
the term W, provided all of them are equal. For in- 
stance, de E G z (p™°) = degG z , degG^p™ ® 

= ( 2 #)( n 2 T)E™=o(r i )(rf-^)( last 
formula is only true for k = Z, otherwise elements con- 
stituting Gz{W) are not all equal). Restricting to the 

case k = Z and denoting W m — P^f'^ <8> (p±) > 
have 



r(w m ) = r. 



/n—m\ 
V 2Z ) 

Gz) ' 



(34) 



It is worth to note that the element gz C Gz is simply 
given by 



,9z=^o[ J56 
n I \n ' 



F 



(n,k) 



2r 



f \g y 



1 7 V n - 1 



(n-l,fc) 



2r 



n-2 



(n-2,k) 



92 



(35) 



(n,fe) _ (m-1) 



M-fc) ( y-i-i \ ( q-p. \ 



with q = 1 — p and gfj 

is the element of a p^* 1 ' corresponding to a block with 
A: particles (the factors r m are due to the averaging while 

the factors come from I pi J ). Restricting to the 
case k = Z, and taking into account 



( F ) 

\n—m/ 



n\ (1-2/Li)' 



(£) (n-m)! F" 
so that 

(n-J (2r\ J_ r 
( L ) \m) 2 m ^ 



(2a*)' 



■L T ' 



(36) 



Ar(i-2 M )™- 



we finally obtain, using (|35p . that 



i-2Z 



m=0 ^ ' 



n-2Z 



= (p-p) n - Z (q-p) Z £ 



p m fn - 2Z 
(p — A*) m V 171 



(p-p) n - z ( q -py 



p 



p- a* 

z 



i-2Z 



(37) 



with go the diagonal element in the same block k = Z . 
In the last calculation we used the relation ^"'^„ 2 f - — 

\n-Z) 

(™ _2Z ) . This proves formula (flU)) for the particular case 
k = Z and arbitrary n. 

For arbitrary k, one proceeds in similar manner as 
for the case k = Z case. Since the respective calculations 
are tedious and not particularly illuminating, we omit 
them and give only the final result: 

n — 2Z min(ro,fe— Z) 

9z = £ M m £ (p~p) n - m - k+l 

m=0 i=max(0,Z+f:-n+m) 

<«-">" C « z ) ("»-")• < 38) 

which, after some algebraic manipulation, can be rewrit- 
ten in the form 



n-k-Z 



gz 



{p-V) n - k {q~p) k £ 



k-Z 

E 



A* 



q-fi 



3=0 

k-Z 
i 



p- a* 



n — k — Z 



= rfp n ~ k q k . 



(39) 



This concludes the proof of the theorem in the thermo- 
dynamic limit L — > oo. For finite L the expression of 
(I18p can be obtained by solving a recurrence relation ob- 
tained from the analytical expressions obtained from (|30|) 
for the cases n = 2,3,4. The correctness of this expres- 
sion can be checked indirectly through the expression for 
the eigenvalues (see [lljV 



VI. SPECTRAL PROPERTIES OF RDM 

Since RDM is block diagonalized with respect to quan- 
tum number k into blocks Bk of size (?) x (?) , to di- 
agonalize the complete RDM we can diagonalize each 
block Bk, k = 0,1, ...,n separately. As mentioned be- 
fore, the permutational invariance of the subsystem of 
size n permits to further block diagonalize each block Bk 
with respect to a quantum number s which labels the 
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irreducible representations of S n (e.g. the YT of type 
{n — s, s}) which are compatible with the block polar- 
ization. This implies that for a block fc the number s 
can take only the values s — 0, ...min(fc,7i — fc). We re- 
fer to the respective eigenvalues as A( s )(ro, fc), omitting 
for brevity the explicit dependence on L,N,r. Note that 
the eigenvalues A( s )(n, k) have degeneracies (™) — ( 
which are equal to the dimension of the corresponding 
S n representation (e.g. the dimension of the YT of type 
{n—s,s}) [Til] - These properties can be easily checked on 
the particular example given in Fig. [5]). By diagonaliz- 
ing matrices Bk for small n one finds that all eigenvalues 
A( s j (n, k) are linear combinations of the matrix elements 
gz of the form: 



A( s )(n,fc) 



min(fc,n — k) 

E ■ 

z=o 



(n, k)g 2 



(40) 



with ccg , integer coefficients, satisfying the following 
properties: 



a s) (n, fc) = 1, 
min(fc,n— k) 



E a { z\n,k) = 



(41) 
(42) 



(43) 



min(fc,n— fe) 

a ( z\n 1 k) = Q fors>0, 
z=o 

min(fe,n— k) f Z\ 

J2 4"' («.*)( ) =0, p = 0,1, ...s- 1,(44) 
z=o \P ' 

min(fc,n — k) 



aW(n,fc)l" J = 

z=o 



(-1)' 



n — k — s 



k — s 



(45) 
(46) 



Note that the coefficients a do not depend on the char- 
acteristics of the original state L, N and r. The depen- 
dence of the RDM eigenvalues on these parameters en- 
ters through the elements gz (fT5|) . This implies that the 
above properties of a can be proved using special cases. 
E.g. in case when all gz = go ; we & r e back to ground 
state problem solved in [l(|. Each block Bk has single 
nonzero nondegenerate eigenvalue Xmj(n, k) — (f£)go, en- 
tailing @3),@g. 

The property (|41 [) is easily proved in the thermody- 
namic limit. In this limit, we have n — 0. Then, all 
gz = 0, for Z > and all -B/c are diagonal implying 
A( s )(n, fc) = .go for any s. This gives fllTl) . 

We have no proof for the properties (f44|) . (|46 [) . Note 
that the (|43|) is a particular case of (|44|) for p = 0. 

Given the exact form of the coefficients or^ (n, k) for 
Z = Q,k, (|4ip and (|4"6")) . the remaining coefficients are 
obtained recursively from (|4^|) . E.g., the last eigenvalue 
of block Bk with the degeneracy (?) — (,,",) is obtained 



using recursively (]44[) for s = k, and p = s 
which gives 



l,s-2,..0 



A (fc) (n,fc) = ^(-l) z [f) 3z , 4 fc) (n, 



k) = {-iy 



(47) 

Note that the last eigenvalue does not depend on n, but 
only on k. Other eigenvalues do depend on n as well. For 
the last but one eigenvalue we have 



4 fc - 1 )(n,fc) = (-l)^(( 



fc-l\ fn-2k + l\fk-l 

Z J { 1 )\Z-1 

(48) 

for Z = 0, 1, ..k. Analogously one can check that 



4*- a) (n,*) = (-l) a [(*- 2 ) 



fc- 2 
Z 

n-2fc + 2\/fc-2 
1 J U - 1 
fc-3 
Z 

n-2k + 3\ fk-3 

Z -2 



4 fc " 3) (n,fc) = (-l)^ 



+3 



n-2fc + 2\/fc-2 
2 J U - 2. 
n-2fc + 3\ /fc-3 

i JU-i 

n-2fc + 3\ /fc-3 
3 



From these expressions the existence of a recursive rela- 
tion is evident, this leading to the following general result 



(n,fc) = (-l) z £(-l) 



m\ [ n — 2k + m 



(k—m) 



with m = 0, 1, fc. Using (09]) together with (gOJ), (fl9j) . 
<|18[) , one obtains the complete spectrum of the RDM. 



k — m 
Z -i 
(49) 



\k-l/ 



VII. CONCLUSION 

In summary, we have provided explicit analytical ex- 
pression of the reduced density matrix of a subsystem 
of arbitrary size n of a permutational invariant quantum 
many body system of arbitrary size L and characterized 
by a state of arbitrary permutational symmetry. We have 
shown, on the specific example of the spin 1/2 Heisen- 
berg model, that the RDM acquires a block diagonal form 
with respect to the quantum number fc fixing the polar- 
ization in the subsystem (conservation of S z or conser- 
vation of the number of particles for non spin models ) 
and with respect to the irreducible representations of the 
S n group. The main results of the paper are represented 
by Eqs. (HHJ), (19]) and ([49]) given above. These results, 
being based only on the permutational invariance and on 
the conservation of S z (number of particles for non spin 
systems), should apply to all quantum many-body sys- 
tems which possess these symmetries and for which the 
mean field theory becomes exact. A more detailed analy- 
sis of the properties of RDM spectrum and of its physical 
implications will be given elsewhere fl2j . 
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